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We discuss results obtained in the study of shape coexistence in the neutron-deficient 72_78 Kr 
isotopes. The method that we used is based on the mixing of axial mean-field configurations after 
their projection on particle-number and angular-momentum. The calculations are performed with a 
Skyrme interaction and a density-dependent pairing interaction. While our calculation reproduces 
qualitatively and quantitatively many of the global features of these nuclei, such as coexistence of 
prolate and oblate shapes, their strong mixing at low angular momentum, and the deformation of 
collective bands, the ordering of our calculated low-lying levels is at variance with experiment. We 
analyse the role of the single-particle spectrum of the underlying mean-field for the spectrum of 
collective excitations. 

PACS numbers: 21.60.-n 21.60.Jz, 21.10.-k, 21.10.Re 21.10.Ky 27.50. +e 



I. INTRODUCTION 

In a mean- field framework , the presence of low- lying 
+ states in the spectrum of an even-even nucleus is usu- 
ally interpreted as the manifestation of shape coexistence 
@: each of the 0+ states, including the ground state, 
corresponds to a mean-field configuration of a different 
shape. Nuclei around the neutron deficient Kr and Sr 
have very early been considered as among the most fa- 
vorable ones for the presence of shape coexistence. First 
studies were performed for Kr and a few neighbor- 
ing odd nuclei with the help of the Nilsson-Strutinsky 
approach Detailed calculations have been carried 
out since then with an improved microscopic-macroscopic 
model and with self-consistent mean-field models us- 
ing non-relativistic Skyrme [f| and Gogny [(| interac- 
tions, as well as relativistic Lagrangians p|- They con- 
firm the presence of oblate and prolate minima in the 
deformation energy surface of some light Kr, Sr and Zr 
isotopes. 

The existence of nearly-degenerate structures corre- 
sponding to different deformations raises immediately the 
question of their stability against dynamical effets be- 
yond a mean-field approach. One can expect that the 
physical states result in fact from a mixing of states with 
spherical, prolate and oblate deformations. Such mixings 
were obtained from models with parameters specifically 
adjusted to the data: the proton-neutron Interacting Bo- 
son Model (IBA-2) and a Bohr-Hamiltonian calcula- 
tion built on a microscopic-macroscopic model 9] . 

An alternative to a mean-field description of shape co- 
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existence is given by the shell model. However, the num- 
ber of active particles and holes necessary to describe 
the neutron-deficient Kr isotopes is prohibitively large, 
and this mass region is out of reach of the conventional 
shell model. The problem is tractable in models per- 
forming a clever truncation of the configuration space 
as in the complex excited VAMPIR approach using a 
phenomenologically modified nuclear matter G matrix as 
residual interaction [ToL EH , or in the Shell-Model Monte 
Carlo (SMMC) method. A SMMC calculation using a 
schematic pairing+quadrupolc interaction is presented in 
Ref. [12. 

First experimental evidence for ground-state defor- 
mation in neutron-deficient Kr isotopes was found in 
Ref. |13|. Subsequent experiments revealed the rich 
structure of the low-lying excitation spectrum in these 
nuclei, with coexisting and mixed bands in 72_78 Kr 
GIGS El. Many more data on transition probabili- 
ties, both in and out of bands , h ave also now become 
available [H El El HI El El IH El El E3- 

The description of such nuclei is a challenge for nu- 
clear structure models. In general, their excitation spec- 
trum varies very rapidly along a given isotopic line. 
For this reason, they constitute a testing ground for 
beyond-mean-field approaches built on Skyrme-Hartree- 
Fock E3- I n this paper, we present an application 
of an angular momentum projected generator coordinate 
method to the description of the low-energy excitations in 
72_78 Kr. At present, our implementation of the method 
is limited to axial mean-field states which are invariant 
under spatial inversion and time reversal. Applications 
of the same method to shape coexistence in neutron- 
deficient Pb isotopes were presented in Refs. [30l l3l| . A 
good qualitative agreement with the data was obtained, 
in particular for the relative position of the coexisting 
bands which result from a mixing of oblate, spherical and 
prolate configurations. However, the mixing is extremely 
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sensitive to many details of the model and on the effective 
interaction, preventing a detailed quantitative agreement 
with the data. Compared to the neutron-deficient Pb iso- 
topes, the description of the neutron-deficient Kr isotopes 
looks simpler. The coexisting bands result from the mix- 
ing of only two kinds of structures, oblate and prolate. 
The number of spherical j shells which contribute to the 
shell structure is smaller than for Pb, and the shells are 
the same for neutrons and protons. It is therefore easier 
to relate the collective states to the underlying mean- 
field. 

Our calculated values presented below have already 
been used for comparison in the experimental report of 
Refs. El El. 



II. THE METHOD 

The starting point of our method is a set of axial 
HF+BCS wave functions. They are generated by self- 
consistent mean-field calculations on a three-dimensional 
mesh in coordinate space 0, l32| . with a constraint on 
a collective coordinate, the axial quadrupole moment 
(Q2o)- I n a spherical nuclear shell model picture, such 
mean-field states incorporate particle-particle correla- 
tions through pairing and many-particle many-hole cor- 
relations through nuclear deformation. As a result, the 
mean-field states break several symmetries of the ex- 
act many-body states. These symmetry violations make 
difficult the connection between mean-field results, ex- 
pressed in the intrinsic frame of reference of the nu- 
cleus, and spectroscopic data, obtained in the labora- 
tory frame of reference. This motivates the second step 
of our method, the restoration of the symmetries associ- 
ated with particle numbers and rotation. Another am- 
biguity in the interpretation of mean-field results arises 
when the deformation energy varies slowly as a function 
of a shape degree of freedom, in particular when a po- 
tential energy surface presents several minima separated 
by a low barrier. In such a case, to assign a physical 
state to each minimum is not a well-founded approxima- 
tion. This problem is eliminated by the third step of our 
method: for each angular momentum, states obtained by 
projecting mean-field configurations of different deforma- 
tion are used as the generating functions of the genera- 
tor coordinate method (GCM). The weight coefficients 
of the mixing are determined by varying the energy and 
by solving the Hill- Wheeler- Griffin equation [3j|. The 
configuration mixing removes the contribution of vibra- 
tional excitations from the ground state wave function. 
It permits also to construct a spectrum of excited states. 
A detailed introduction to the method has been given in 
Ref 0. 

The same effective interaction is used to generate the 
mean-field states and to perform the configuration mix- 
ing: the Skyrme interaction SLy6 |36| in the particle- 
hole channel and a density-dependent, zero-range inter- 
action in the pairing channel. As required by the SLy6 



parametrization, the full two-body center-of-mass correc- 
tion is included in the variational equations to generate 
the mean-field and in the calculation of the projected 
GCM energies. The strength of the pairing interaction is 
the same as in previous studies of light and medium-mass 
nuclei, i.e. —1000 MeV fm 3 for protons and neutrons in 
connection with cutoffs above and below the Fermi en- 
ergy, as defined in Ref. [37] . Similar implementations of 
angular-momentum projected GCM have been developed 
also for the non-relativistic Gogny interaction [35] and for 
relativistic mean-field Lagrangians [38l| . 

To define a dimensionless deformation parameter, we 
use the sharp edge liquid drop relation between the axial 
quadrupole moment and a parameter 
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with R = 1.2 A 1 / 3 fm. The mass dependence of the 
quadrupole moment is partly removed in this parameter 
/?2 • It should not be compared with the multipole expan- 
sion parameters that are used in microscopic-macroscopic 
models whose origin is not the same and usually are sig- 
nificantly smaller. 

With our method, B(E2, J — > J') values for in-band 
and out-of-band transitions as well as spectroscopic mul- 
tipole moments Q S (J) are determined directly in the lab- 
oratory frame of reference [3^ . As we use the full model 
space of occupied states, we neither have to distinguish 
between valence particles and a core, nor is there a need 
for effective charges. 

To make a comparison with other approaches easier, 
it is useful to define quantities similar to intrinsic frame 
deformations from the spectroscopic and transition mo- 
ments. An intrinsic transition charge quadrupole mo- 
ment Q^' can be derived from the B(E2) values within 
the static rotor model 
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In the same way, an intrinsic charge quadrupole moment 
can be related to the spectroscopic quadrupole moment 
Q C (J) via the relation 

For axially-symmetric shapes, i.e. pure K — states, this 
relation reduces to Q {s) {J) = ~Q C {J) (2J + 3)/J, with 
a change of sign between and Q C {J). Dimensionless 
deformation parameters f3^ and (3^ can be determined 
from Qw anc j Q( s \ respectively, through an equation 
similar to Eq. Q), with A replaced by Z. Within a given 
band, and might still depend strongly on angu- 
lar momentum. 

Q( s ) and measure different properties: de- 
pends on a single state, while Qw probes the geometry of 



3 




-0.4 -0.2 0.0 0.2 0.4 0.6 



02 

FIG. 1: Nilsson plot of neutron single-particle energies with 
positive (solid lines) and negative (dotted lines) parity as a 
function of the intrinsic mass deformation fxf \ as obtained 
for 74 Kr. 



the initial and final states whose wave functions can cor- 
respond to very different mixings of projected mean-field 
states. In general, these two moments take different val- 
ues. Their near equality indicates that the assumptions 
of the static-rotor model are fulfilled, i.e. a well-deformed 
rotational band not mixed with other bands. 



III. RESULTS 
A. Potential landscapes 

Figure ^ provides the Nilsson diagram of the self- 
consistent single-particle energies for neutrons. The di- 
agram for protons looks very similar, but shifted in en- 
ergy. One can see two main differences with the single- 
particle energies used in some other methods. In the 
Woods-Saxon and the Folded- Yukawa potentials used in 
the microscopic-macroscopic calculations of Ref. Q and 
[40| respectively, as in the Woods-Saxon potential used 
to generate the single-particle spectrum for the SMMC 
calculations of Ref. [12||. thep 3 / 2 andpi/ 2 levels are closer 
to the <7g/ 2 level by approximately 2 MeV and the sepa- 
ration between the /s/ 2 and the <?g/ 2 levels is larger. As 
a result, the p 3 / 2 and the / 5 / 2 levels are much closer and 
their order is even changed in the case of the latter two 
models. These modifications have a strong effect on the 
deformed gaps which may correspond to quite different 
deformations and vary in size. The spin-orbit splitting 
for the / and p levels in the three potentials are very 
similar to ours. The differences between the single parti- 
cle spectra must thus be related to the relative position 
of states with different orbital angular momentum within 
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FIG. 2: Mean-field (dotted) and J = projected deformation 
energy curves (solid) for r2_78 Kr (see text). 

a given shell. 

Figure [3 shows the mean-field and J = projected en- 
ergy curves for 72 ~ 78 Kr obtained with the SLy6 Skyrme 

parametrization. Throughout this paper, the projected 

(i) 

energy is plotted as a function of the intrinsic /3 2 value 
of the mean-field states from which they are obtained. In 
our opinion, this quantity provides the most convenient 
and intuitive label that can be defined for all states, irre- 
spective on the level of modelling. However, it should not 
be misinterpreted as an observable. With our method, 
one calculates transition and spectroscopic multipole mo- 
ments in the laboratory frame, which can be directly 
compared to experimental data. However, the spectro- 
scopic moments which characterize a state do not pro- 
vide useful coordinates to plot potential energy curves as 
they scale with angular momentum and are even identi- 
cally zero for J — 0. For projected states with J > and 
large intrinsic deformation in nuclei with A larger than 
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100, 02 i s very close to the intrinsic deformation (3^ 
determined from the laboratory frame quadrupole mo- 
ment through the static rotor model, Eq. Note 
that any coordinate might be misleading for the GCM, as 
the metric is related to the inverse of the overlap matrix, 
which has no direct connection to any deformation. 

The mean-field energy landscapes (dotted lines in Fig. 
|2J) show that the energies of the four nuclei vary very 
slowly with deformation. Our calculations predict also 
a transition from a nucleus with coexisting prolate and 
oblate minima in 72 Kr to a soft, spherical, anharmonic 
vibrator in 78 Kr. The many shallow local minima and 
plateaus in the total energy curves can be directly related 
to the gaps in the Nilsson diagram in Fig. 

The two minima in the mean-field energy curve (dotted 
line) of 74 Kr reflect the N — 38 gaps in the Nilsson dia- 
gram at small oblate and large prolate deformations. For 
76 Kr, the spherical mean-field minimum is related to the 
large spherical N = 40 subshell closure, while the shallow 
oblate and prolate structures correspond to two deformed 
N = 40 gaps in the Nilsson diagram. The prolate mini- 
mum has moved to smaller deformation than in 74 Kr, the 
deformed N — 38 gap being at larger deformations than 
the N = 40 gap. The spherical N = 40 subshell closure 
is strong enough to stabilize the spherical shape up to 
the N = 42 isotope 78 Kr. 

As can be seen from the solid lines in Fig. El the energy 
landscapes are qualitatively modified by the projection 
on J = 0. Since a spherical mean-field state is already 
a J = state, the energy gain by projecting it on angu- 
lar momentum is zero while the projection of a deformed 
mean-field state always leads to an energy gain, which 
increases very rapidly at small deformation. In almost 
all spherical and soft nuclei, this creates minima at pro- 
late and oblate deformations with j/J^ | values around 

0. 1. These states usually have a large overlap close to 

1, which means they are not two different states, but 
represent the same one, which can be associated with a 
"correlated spherical state" . In nuclei with shallow mean- 
field minima at small deformations as the Kr isotopes, 
the projection merges this near-spherical spherical mini- 
mum with the slightly oblate one into a broad structure. 
The prolate mean-field minimum being at larger defor- 
mation, two distinct structures appear on the prolate side 
in 74 ~ 78 Kr, a "spherical" one at a value around 0.15, 
and a well-deformed one at a @2 value around 0.5. 

The topology of the 72 Kr energy surface on the oblate 
side can hardly be explained by a single large oblate 
N = 36 gap in the Nilsson diagram. The mean-field and 
J = minima at large oblate deformations seem to result 
from a delicate balance of shell effects from N = 34 to 
N — 40 which reduce the level density around the Fermi 
surface sufficiently to create a shallow oblate minimum 
and a very soft energy surface. The shallow minimum on 
the prolate side can be associated with two close N = 36 
gaps. 

The fact that the structures in the potential land- 
scapes can be easily associated with gaps in a Nilsson 



diagram is an attractive feature of our approach which 
makes straightforward the connection with simpler mod- 
els. This is an advantage compared to the interact- 
ing shell model, in which states are constructed as non- 
intuitive np-nh states in a spherical basis. 

One can wonder whether it is meaningful to restrict the 
description of the light Kr isotopes to axial shapes. There 
have been a few explorations of the triaxial degree of free- 
dom which can help to answer this question. Yamagami 
et al. ^3 na ve calculated the potential landscape of 72 Kr 
along several shape degrees of freedom, including non- 
axial octupole deformations, using the Skyrme SIII inter- 
action. Their potential energy curve is similar to ours for 
purely quadrupole axial deformations. They found that 
the prolate and oblate energy minima are separated by 
a triaxial barrier of 500 keV, and that the oblate shapes 
are soft with respect to non-axial octupole deformations. 
Bonche et al. have obtained rather large triaxial barri- 
ers for 74 Kr and 76 Kr, also using the SIII interaction. Let 
us also mention that Almehed and Walet ^3 have self- 
consistently determined a collective path in 72 Kr for dif- 
ferent values of the angular momentum in a small model 
space using a schematic interaction. For J = 0, they ob- 
tain a purely axial path that connects the oblate ground 
state and the prolate minimum. For J = 2, the path 
that they obtain does not cross the spherical configura- 
tion and makes a small excursion through triaxiality, as 
can be expected from the projected energy curves. From 
these studies of triaxiality, one can conclude that in most 
models, the oblate and prolate minima found as a func- 
tion of the axial quadrupole moment are true minima. 
The barrier between the minima in the 7 degree of free- 
dom may not be very high, but to restrict the model to 
axial symmetry can be expected to provide a valid first 
approximation. 



B. Spectroscopy 

We first focus on the three isotopes for which very 
complete sets of experimental data are available. This 
enables us to introduce our notations and all the ingre- 
dients which are necessary for a comparison between our 
results and the data. We will come back to the more 
prospective case of 72 Kr at the end. 



In Fig. are given the lowest projected GCM levels 
for J = up to 8 together with the corresponding pro- 
jected energy curves. The collective levels are plotted at 
the mean deformation of the mean-field states from 
which they are constructed, defined as 

4 l) = E^ 9lM% (4) 
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FIG. 3: (color online). Mean-field energy curve, angular- 
momentum projected energy curves and projected GCM for 
74 Kr. The projected energy curves are plotted as a function of 
the intrinsic deformation of the mean-field they are projected 
from, the projected GCM levels at the average deformation 
/? 2 , Eq. (QJ. The labels on the right-hand side give the total 
binding energy as calculated, the labels on the left-hand side 
the energy relative to the J = GCM ground state. 



which provides in many cases an intuitive picture of the 
band structure in a nucleus. Figure however, illus- 
trates the limits of the meaning of /?2- In this figure are 
shown the collective wave functions gf for J = 0, 2 and 
4. All low-lying + states result from a strong mixing be- 
tween prolate and oblate mean- field configurations. The 
values of /3 2 for the ground state and the second excited 
state reflect the dominance of oblate and prolate defor- 
mations. For the ground state + state, the very small 
value of 02 does not mean that this state is nearly spheri- 
cal, but that the weights of prolate and oblate shapes are 
nearly equal. For higher J values, the mixing between 
oblate and prolate configurations is less pronounced and 
the value of fa represents better the structure of the 
states. 

As could be seen in Fig. |3 the J = projected energy 
curve for 74 Kr presents two nearly degenerate minima at 
oblate and large prolate deformations, and a third shal- 
low excited minimum at small deformation. This struc- 
ture is reflected in the spectrum of GCM states plotted 
in Fig. [3] and in the wave functions drawn in Fig. 0] 

As can be seen in Fig. 21 the J = GCM wave func- 
tions strongly mix projected oblate, near-spherical and 
prolate configurations. The near degeneracy of the oblate 
and prolate minima is lifted, the first excited + state 
having an excitation energy of 0.49 MeV. The situation 
becomes simpler for larger angular momenta. The first 
2 + level still strongly mixes prolate and oblate deforma- 
tions and has a wave function rather similar to the first 
excited + state. On the contrary, the second 2 + state 
is mainly constructed from oblate mean-field configura- 
tions projected on J = 2 with a negative mean deforma- 
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FIG. 4: (color online). Collective wave functions (?j7' for 
the lowest states with J = 0, 2, 4, and 6 for 74 Kr as a func- 
tion of the intrinsic deformation of the mean-field states they 
are constructed from. The corresponding energies and aver- 
age deformations 02 are plotted together with the projected 
energy surface. 



tion larger in modulus than that of the ground state. The 
two first 4 + and 6 + states are dominated by only prolate 
or oblate deformations. One can thus expect that an ap- 
proximation based on a fixed mean-field configuration, 
like the cranked mean-field methods, is justified beyond 
J = 4. 

The nodal structure of the collective wave function is 
more complicated than the structure that one would ob- 
tain in a one-dimensional potential well. It reflects the 
fact that the rotation of a deformed wave function around 
the j/-axis generates wave functions which cannot be rep- 
resented along a single axis: the exchange of the x and z 
axes generates prolate and oblate shapes for two differ- 
ent values of the angle 7. For this reason, the angular- 
momentum projected GCM wave functions is not one- 
dimensional, even when only axial states are mixed |39| . 

74 Kr provides an excellent example that a compari- 
son between the intrinsic deformation of a mean-field 
minimum and the transition quadrupole moment derived 
from a B(E2) value might be misleading. There are three 
reasons for that: (i) for light and medium-mass nuclei, 
the minima of the projected energy curve are usually at 
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FIG. 5: Spectrum of collective states for 74 Kr as seen in 
experiment (left) and as obtained from our calculation (right). 
The arrows with numbers denote the B(E2, J — > J — 2) values 
for transitions between states with different J and B(E2, J — > 
J) for transitions between states with the same J, both given 
in e 2 fm 4 . Experimental values for the B(E2) values are taken 
from IHEIIl. Tne labels "oblate" and "prolate" given to the 
bands correspond to the main components of the collective 
wave functions. 



larger prolate and oblate deformations than the minima 
of the mean-field energy curve, (ii) the minima of the 
projected J = and J = 2 energy curves correspond to 
the projection of mean-field states with different defor- 
mation and the assumption of a static rotor is not valid 
at low angular momentum, and (iii) the wave functions 
of the low-lying states are quite broad and strongly mix 
oblate or prolate configurations; moreover, the mixing is 
different for the + states and 2 + states. 

A schematic two-level-mixing model has been used by 
Becker et al. [TEl |43| to analyze the experimental data. 
In this model, the energies of "pure" prolate and oblate 
configurations are deduced from an extrapolation to low 
spins of the high spin parts of the bands. This is valid 
if configuration mixing decreases with spin, which seems 
reasonable in view of our results. Data are described well 
by assuming that the unperturbed levels are nearly de- 
generate with only a 20 keV energy difference. This leads 
to nearly equal contributions from the oblate and prolate 
unperturbed states. This result is consistent with our cal- 
culations, although the assumption that one has to mix 
only two configurations is too crude. The complicated 
structures that we obtain for the collective wave func- 
tions, of Fig. 0] cannot be reduced to the mixing of two 
or even three configurations with a well defined shape. It 
does not seem therefore possible to connect our results 
with a phenomenological two-level mixing model. 

Theoretical and experimental excitation spectra and 
B(E2) values are compared in Fig. [5] As in the experi- 
ment, the theoretical levels are ordered in bands on the 
basis of their spectroscopic quadrupole moments and of 
the dominant E2 transitions. A low- lying 7 band, which 
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FIG. 6: (color online). Comparison of calculated and ex- 
perimental Q^ 3 \j) and Q^\j — > J — 2) moments for states 
in the prolate and oblate bands, assuming that all states are 
purely axial. Note that —Q^' is plotted for oblate states for 
better comparison with Q^ s \ Experimental values are taken 
from El. 



has a 2 + band head at 1.74 MeV, has been omitted since 
K = 2 states are outside our model space. The B(E2) 
values are always given for transitions from J — > J — 2 or 
from J — > J, although values found in the literature cor- 
respond sometimes to B(E2, J — » J+2) values, which are 
by definition larger by a factor 5. Our calculation is only 
partly consistent with experiment. We reproduce the co- 
existence of oblate and prolate bands, strongly mixed 
at low angular momentum, but the lowest + states do 
not have the right ordering. Our calculated ground state 
is predominantly oblate, but the prolate band becomes 
yrast already at J = 2. For higher spins the spectrum 
is too stretched out, a problem that has already been 
encountered in previous GCM studies 0, ^ . 

The in-band transition probabilities are slightly over- 
estimated for both bands, which hints at a slightly 
too large deformation of the dominating configurations 
in both bands. Nevertheless, all calculated values are 
within a factor two of the experimental ones. Out-of- 
band transitions between states at the bottom of the 
bands are large, in particular the experimental value of 



B(E2,2+ IO 0+ bl ) of 938+^ u e 2 fm 4 given in |25j, or 
1120 ± 460 e 2 fm 4 given in 16]. Our predicted out-of- 
band B(E2) values become weak already at J = 4. 

The theoretical spectroscopic quadrupole moments are 
compared with experiment in Fig.|SJ These moments pro- 
vide by their sign the only unambiguous way to assign 
a prolate or an oblate character to a band. To com- 
pare them with in-band B(E2) values as well, and 
Q(') values derived from the static rotor model assuming 
K = are shown. Except for the oblate 2 + state, our 
calculated and Q^' values are very similar. Exper- 
imental values for and Q^' are also similar within 
the large error bars of Q^ s \ Our calculation gives sys- 
tematically larger deformation for both the prolate and 
oblate states as already discussed above. 
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2. 7b Kr 

The projected energy curves in 76 Kr, plotted in Fig. [7] 
are very similar to those obtained for 74 Kr, with coex- 
isting prolate and oblate bands, except that the prolate 
minima at fa values around 0.45 are less pronounced. As 
a consequence, the ground state wave function is more 



dominated by oblate configurations than that of 74 Kr, as 
confirmed by the shape of the collective wave functions 
given in Fig. This is at variance with experiment, 
which indicates a prolate ground state, see also Fig. [S] 
Again, we omit from Fig. the experimental states as- 
signed to a 7 band with a 2 + band head at 1.222 MeV; 
see, e.g. 0,H^]. The calculated prolate band head has 
a higher excitation energy than for 74 Kr. The prolate 
band becomes yrast at J = 4. Compared to 74 Kr, the in- 
creased "purity" of states within the prolate and oblate 
bands happens at higher angular momenta, which is re- 
flected in the large out-of-band B(E2) values up to J = 4 
in Fig. El 

Figure ^| compares the intrinsic quadrupole moments. 
As for 74 Kr, our calculated values for and Q^> are 
very similar, with small deviations at low J. Notably, this 
is not the case for the experimental values when assuming 
K = 0, which might hint at a small admixture of triaxial 
shapes to the prolate band, which is outside the scope of 
our calculation. Our calculated Q^> are very close to the 




FIG. 10: (color online). Same as Fig.© but for 76 Kr. Ex- 
perimental values are taken from |25ll . 
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FIG. 11: (color online). Same as Fig.EI but for 78 Kr. 
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FIG. 13: The same as Fig. but for 78 Kr. Experimental 
values for the B(E2) values are taken from [2(J. 
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FIG. 12: (color online). Same as Fig.EJ but for 78 Kr. 



experimental values for the prolate band, while the 
of the 2 + state in the oblate band is underestimated. 



Among the Kr isotopes, 78 Kr has the lowest excitation 
energy for the first 2 + state, and the largest B(E2, 2^ — » 
0i) value, corresponding to a transition quadrupole de- 
formation (3^ of 0.35. The rotational band built on the 
ground state has been seen up to J = 24. Two additional 
bands are known at low excitation energy, a first usually 
interpreted as a 7 band built on a low-lying 2 + state 
at 1.148 MeV, and a second one built on a + state at 
1.017 MeV, originally assumed to be oblate. According to 
the recent analysis of [2(|, the spectroscopic quadrupole 
moment of the 2 + state in this band is negative, which 
indicates that this band should, in fact, be prolate. 

Looking to the energy curves of 78 Kr, the collective 
states obtained from the projected GCM, Fig. ^2 an d 
the collective wave functions, Fig.^J one can attribute 
to the ground state a dominating oblate structure, and to 
the band head of the first excited band, located at 2 MeV, 
a well-deformed prolate shape. All low-lying states are 
strongly mixed, which leads to large out-of-band B(E2) 
values up to J = 6. 

Data for transition moments are taken from a recent 
Coulomb excitation experiment [2(J. Values for in-band 
transitions in the yrast band were published earlier in 
Refs. |M E3, with a B(E2; 2+ -> 0+) value of 91(5) e 2 
fm 4 in [l7|, compatible with the value of 130(40) obtained 
in 

As can be seen in Figs.ll3land ll4l our values are close 
to the experimental data for the low-lying states in the 
prolate band. The slight increase of with J that 
we obtain can be related to the gradual disappearance of 
the prolate minimum at P2 values around 0.45 and to the 
deepening of the shoulder at 02 values around 0.6, which 
becomes the prolate minimum at J > 8, see Fig. ^2 hi 
contrast, the experimental Qj^'s are constant or slightly 
decreasing with J, see also |2fll27^ . 

The large mixing between the prolate + state with an- 
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FIG. 14: (color online). Same as Fig. | 
perimental values are taken from |2pL 126 



but for 78 Kr. Ex- 



other excited + state considerably reduces the in-band 
B(E2; 2+ ro -> 0+ ro ) value, as the E2 strength from the 
2p ro level is nearly equally distributed. Therefore, our 
B(E2; 2+ ro — » 0p ro ) probability is slightly smaller than 
the experimental B(E2; 2f — > Oj) from the correspond- 
ing experimental state. As a consequence of this strong 
mixing, the calculated and Q" values differ more 
than for 74 Kr and 76 Kr, QW being larger than Q^> ex- 
cept for J = 4. Experimentally, values are smaller 
than QW ones 0. 



4- 



2 & 



Rotational bands at high spin in 72 Kr have been in- 
vestigated extensively, mainly motivated by the quest 
for fingerprints of T — proton- neutron pairing [44l l45| . 
Much less is known about the structure of this nucleus 
at low excitation energy. Only one band and two states 
that cannot be grouped into bands have been observed 
at low spin 0|. The precise structure of the low-spin 
yrast states is not completely clear either, as the ground 
state is argued to be oblate, in parts by consistency with 
theoretical predictions, while at higher spin the yrast 
states are assumed to be prolate. It was also argued 
in Ref. that the large value obtained from the 
B(E2;8+ -> 6+) = 2090(780) e 2 fm 4 favors its inter- 
pretation as a transition between prolate states. 1 We 
obtain, however, a much larger B(E2) for the 8p ro — ► 6+ ro 
transition than the experimental value, which is in fact 
consistent with our 8^ bl — > 6^ value. In view of the sys- 
tematic overestimation of the experimental B(E2)s at 
large spins given by our method for the heavier Kr iso- 
topes, one cannot draw a conclusion on the nature of the 



1 The authors of Ref. use a definition of fii different from ours, 
so that the values cannot and should not be compared. 
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FIG. 15: (color online). Same as Fig.[2for 72 Kr. 
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FIG. 16: (color online). The same as Fig. H for 72 Kr. 



8 + state in 72 Kr. However, it is clear that these large 
B(E2) values do not exclude an oblate band as previ- 
ously argued. 

The potential energy curves, together with the GCM 
spectra, are shown in Fig. El and the collective wave 
functions in Fig. 1161 According to the usual interpre- 
tation of the low-lying + states, our calculation repro- 
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FIG. 17: The same as Fie. ©for 72 Kr. There is no unambigu- 
ous assignment of the experimental yrast states into a prolate 
or oblate band (see text). Experimental data for the B(E2) 
values are taken from (8 + -> 6 + ) and HI (2+ ^ 0+). 



duces their order with a dominantly oblate ground state. 
In Fig. 1171 the experimental and theoretical excitation 
energies and B(E2) values are compared. According to 
the usual interpretation of the two low-lying + states, 
our calculation reproduces their order with a dominantly 
oblate ground state. The first excited 2 + is at about the 
same excitation energy as the experimental one, while 
the first excited + is slightly too high. Our calculated 
oblate states are yrast up to J = 6. With the excep- 
tion the low-lying + states, which strongly mix oblate, 
spherical and prolate states, the mixing between the pro- 
late and oblate configurations is less important than for 
the heavier Kr isotopes. Within the accuracy that can 
be expected from our model, our B(E2) value for the 
2i, — > 0^~ bl transition is consistent with the experimen- 
tal one, B(E2; 2+ -> 0+) = 1000(130) e 2 fm 4 . 

72 Kr is the only Kr isotope studied here for which we 
can assume to reproduce the ordering of the low-lying 
states. 



C. E0 transitions and radii 

The nuclear matrix element entering the monopole de- 
cay rate of a J 71 ' state to a state with the same spin and 
parity is given by ^(| 



PEo(Jk' — ► Jk) - 



(JMk\f 2 \JMk') 



R 2 



(5) 



where R — 1.2 A 1 / 3 fm. Within a simple model this 
matrix element can be related to the amount of mixing 
of configurations with different deformations in the phys- 
ical states. As it is a non-diagonal matrix element, it is 
also very sensitive to the detailed structure of the initial 
and final states. Within the error bars, the experimen- 
tal p 2 E0 values are very close for 72_76 Kr, and slightly 




36 38 40 42 36 38 40 42 
Neutron Number N 

FIG. 18: (color online). Left panel: comparison of calculated 
values for p|;o(0^~ — > 0^ ) with experiment. Experimental val- 
ues are taken from [3 ( 72 " 74 Kr) and [13 ( 76 - 78 Kr). Right 
panel: isotopic shift of the mean-square charge radii for the 
lowest (dominantly) oblate + state (triangles) and the low- 
est (dominantly) prolate + state (squares) . The ground state 
predicted by our calculation is always oblate. 



smaller for 78 Kr. The variation of the GCM values is 
much larger, but in view of their sensitivity to model de- 
tails it is encouraging that our calculated p E0 values are 
within the experimental error bars for 74 Kr and 78 Kr and 
underestimate the value for 72 Kr by only a factor two. 
The large underestimation for 76 Kr can be related to the 
much smaller mixing between the oblate and prolate + 
states that we obtain for this isotope. 

The reasonable description of the p E0 states indicates 
that, with the exception of 76 Kr, we have about the right 
amount of mixing between the oblate and prolate + 
states, which is independent of the relative placement of 
the energy levels. In contrast, the systematics of ground- 
state radii is not well described, as our calculations give 
the "wrong" ground state for 74_78 Kr. This can be seen 
from the right panel of Fig. ^| where the isotopic shifts 
of the mean-square charge radii 



(6) 



5r 2 (N) = r 2 {N) - r 2 {N = 50), 



are plotted. The Kr radii present an anomaly |48l |49j : 
they are larger for 74 ~ 84 Kr than for the N = 50 isotope 
86 Kr. The standard interpretation is that the ground 
states have an increasing admixture of prolate states of 
large deformation when going to more neutron-deficient 
nuclei, which overcompensates the volume effect of de- 
creasing mass. As our calculations always predict domi- 
nantly oblate ground states, we obtain a negative isotopic 
shift for all for isotopes (triangles in the right panel of 
Fig. I18|) . Our calculated excited dominantly prolate + 
states, by contrast, lead to positive isotopic shifts for all 
four isotopes. It is unclear how well the ground state of 
86 Kr itself is described within our method, so one has to 
be careful with the interpretation of the absolute values 
of the isotopic shifts. However, the Sr 2 data confirm once 
again that the prolate states should have a much larger 
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contribution to the ground state. 

The relativistic mean-field calculations with the NL- 
SH interaction presented in Ref. Q do reproduce the sys- 
tematics of the isotopic shifts, which hints at the more 
realistic single-particle spectrum of this interaction in the 
fp shell. 

IV. DISCUSSION AND CONCLUSIONS 

Our results show that in medium-mass nuclei with co- 
existing shallow minima in soft potential energy land- 
scapes, one has to be very cautious when comparing ex- 
perimental data and the mean-field states corresponding 
to local minima. Projection on good angular momen- 
tum alters significantly the potential energy surfaces. Af- 
ter configuration mixing, the collective wave functions 
present a very large spreading and may extend over sev- 
eral mean- field minima. We also find rotational bands, 
for example in 78 Kr, which are not related to a minimum 
in the mean-field energy landscape, but to a plateau. 

While our calculations reproduce many of the global 
features of the neutron-deficient Kr isotopes, there re- 
main several discrepancies with experiment when looking 
in detail. 

• The collective spectra are too spread out in exci- 
tation energy, which appears to be a general prob- 
lem of GCM of angular-momentum projected axial 
mean- field states | 3 llloCi . This problem could prob- 
ably be cured by a larger variational space for the 
projected GCM, in particular by treating dynam- 
ically the pairing, by breaking axial symmetry or 
time-reversal invariance. Breaking axial symmetry 
will also permit the description of 7 bands. Work 
in these various directions is underway. 

• Our study of neutron-deficient Kr isotopes clearly 
points at a problem with the SLy effective interac- 
tions in the fp shell, related to the relative position 
of a few single-particle levels. 

As already discussed in the literature, (see, e.g., the 
discussion in [l2|), the spectroscopic properties of the 
neutron-deficient Kr region are extremely sensitive to 
the details of the shell structure. The number of lev- 
els around the Fermi energy is quite small which results 
in several shell gaps for every neutron and proton num- 
ber from 34 to 42. Any change in the relative position 
of levels and of their ordering will modify the size and 
deformation of the gaps in Fig. ^ 

Our results, and those of models with adjustable single- 
particle spectra that perform better for Kr isotopes, sug- 
gest a few necessary modifications of the single-particle 
spectra obtained with the SLy6 Skyrme parametrization: 

• The experimental evidence that there are no spher- 
ical + states in the spectrum of 76 Kr shows that 
the N — 40 spherical shell gap is too strong. This 



is consistent with our results for the N ss Z » 40 
region obtained in a systematic study of mass and 
deformation where the ground state for 80 Zr was 
obtained spherical with a very pronounced shell ef- 
fect while experiment suggests that is a well- 
deformed rotor [51|. A g 9 / 2 level closer to the p 1 / 2 
orbital would decrease the N = 40 spherical gap. 
It would also decrease the oblate N = 38 gap, and 
shift the prolate N = 38 gap to smaller deforma- 
tion, probably reducing the deformation of 74 Kr. 

• Taking the Wodds-Saxon single-particle level 
schemes of Fig. 17 in 40] as an example, a de- 
crease of the separation between the and the 
P3/2 levels should change the single-particle spec- 
tra at small prolate deformation: a gap at N = 38 
that extends from oblate shapes to prolate shapes 
would replace the gap at N = 36 around (3 2 — 0.2. 
It can be expected that this would drive 74 Kr to- 
wards prolate shapes. 

We have checked that all modern successful Skyrme inter- 
actions give single-particle spectra in the fp shell which 
are similar to those obtained with SLy6. This points to a 
deficiency of the standard Skyrme interaction in general 
and not to a difficulty intrinsic to SLy6. Other models 
as the relativistic mean-field model (RMF) might seem 
to work more satisfactorily in the fp shell and to pre- 
dict potential landscapes in better qualitative agreement 
with the Kr data than ours 0. However, those models 
do not describe correctly the neutron-deficient Pb region 
|52j,|53| , where many Skyrme forces perform quite well 

Altogether, the present study confirms the conclusions 
drawn in ,50] that the current functional form of the 
Skyrme energy functional is not yet flexible enough to 
cover the relevant aspects of nuclear structure with the 
same good quality for all regions of the nuclear chart. 

The combined analysis of Figs. ^ and [3 suggests that 
the structure of collective states in the neutron-deficient 
Krypton isotopes provides a sensitive testing ground for 
future attempts to construct better effective interactions. 
It has been pointed out recently that a tensor interaction, 
absent in the existing standard mean-field models, in- 
troduces a p article-number dependence of single-particle 
spectra |53, |£g , that might resolve at least some of 
the deficiencies in the fp shell we encounter here, and 
that are known for other regions of the nuclear chart. 
Work in that direction is underway. 
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